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ABSTRACT 

An important example of an event which obeys the geometric 
probability law with parameter 9 is the number of cycles required 
to obtain the first failure of an item, where the success of each 
cycle is independent with probability of success @. 

The true probability of success, in manufactured items, is 
usually unknown and must be estimated on the basis of observed 
data obtained from a test of sample items. An important estimator 
for this probability is the unbiased estimator, defined as that 
estimator whose expected value equals the true value of the para- 
meter. In this study, an unbiased estimator for 9 is derived. 
This estimator is based on the results of a series of independent 
items, each cycled to first failure. 

Series approximations for the variance of this estimator are 
derived, and some values of the variance are tabled for those 


cases thought to be of special interest. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 


Symbol Description 
8 Probability of success of a single cycle 
n Number of items cycled till first failure 
S Total successful cycles completed in n tests 
N, The number of ways n tests can sum to s successes 
8. The unbiased estimator for 9 
M. The product (Ne, 9, ) 
R The ratio =) 
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ine Introduction. 
The geometric probability law with parameter @, where 
0O£6=1, is specified by the probability mass functions: 


p(x) ar (1 - 8) fora = Oni? 


= Q Otherwise, 
where x is the number of successful cycles to first failure of 
one item. 
When n independent items are cycled to first failure, the 


probability of each ordered event is 
a Ax: S 
] Jo{1-0) = @°(1-6)" where s = 2%; 
¢=/ 


The probability of exactly s successes is thus the union of 
all events with exactly s successes. The number of such events is 


the number of ways n tests can sum to s successes. Therefore; 


p(s) 


N.@°(1-6)" for SSS OF 1, Dee Jes esecees 


OQ otherwise. 

Since N, M, and 6 are functions of both n and s, the subscript 
(n,s) is necessary to uniquely identify the value. However, since 
n is a known variable which is fixed prior to testing, the abbre- 


viated subscript s will be used. 


ris ie oi SD AS GG Ne 





Ze Derivation of the Unbiased Estimator. 

By definition, and estimator 6. is unbiased if its expected 
value equals the true value of 0, where the expected value E(@.) 
is defined by 

A ” 
E(@,) = > 6. p(s) over all s such that p(s)>0. 
s 


Thus, 6. is an unbiased estimator for 9 if and only if 
oo 
e-> 6 N. 9° (1-0)" 
a s "Ss (1-0) 
$20 


» M, O° (1-0)" . 


520 


By expanding the (1-0)" term, the summation becomes 


e+ Sy ot [2 - (tos (alot =~ - entered + co? 0] 
: . | 1- (f)e +(S)e? ----- +(-1P7H 2 on"? + (-1)" oF 
= m | 2 - (i)e ‘ (5)e° ee (1h yom + G2)" ett] 
ae 3 (He? + (Z)e" - ~~ -- +1 Bent + CD" had 
+------ ete. 


+ 07] M, - (2x, + ( 3}, | 


. 0] a, - (P)o + (3) - (3) | 


+ === etc., until the @"th term. Then, 





+ 9” ly, - (Ta + - - + ern) a (na) + (1s H, | 


+ oT as = (By to DEE (aaa + CDE] 
ager | M2 SI ygp te = = (-1)9* (2) M3 + CD" Ms | 
. cout coe etc 


By equating coefficients: ES = 0 
n 
I, ; (7) % ] =e 
All remaining coefficients are zero. 
Solving for Ms (s = 0, 1, 2, ---), the following table of equations 


is produced: 


Mo = 0 

M = 1 
- (7 

My = (JM 


. {2 cea ia 
Be i (7) My ( Boy 
Example of table (1) 


My = (2) 05 ; (3), +(3)my with n= 5. 


Ms = (T}my - (2)ms + (3)m> - (8) ! 

Me = (B)Ms - (B)M, + (3)m5 - (2)m + £8) | 

My = (T)Me - (2)Ms + (3)m, ~ (2) + 3) (8) 41 

Me = (El ~ (Be + (S)s5 ~ (2) + (Shea (elm + (Fm 


ot 
+ 


and so forth. 


3 





Initial numerical solutions for M, are as LoLlows < 


M, = 0 

M, = 1 

M, = 

Hs = (So - (B) = aGoetd = (99 

ue (PE) (Ge + (3) = SDD ~ (99) 


a 


Thus , M. can be written as the expression 


= at EOP 60 lta 2 oc 
s-l 


69) 


The proof that this expression holds for all s will be by induction. 
The assumption: M; - Og for 1 = 0, 1, 2, 35> -= oe 


To prove: M, = Ge 


From equations (1), M, can be written as the following summation of 


the M;: 


My 


Scot (lesa) - [Sot (yes? 
[ES cn! (MET) - (EN?) 


i:0 


Mm = — bx be (? He" fe y eal : oo 


Let X = ntk-1, and Y = ntk-2. Then, 


> cot (8) (4) - Sent (2) (54) + (ES?) 
4:0 


(20 


But ) = ca _ head, ‘ Performing this transformation, we have 





By use of the following combinatorial identity: 


S (-1)* (3 (hee) =e fence > TY, 
t= 0 


which is derived in Appendix I, the desired solution is obtained. 


= = n+k-2 = n+k-2 
toe 1 +("2577} aperi 


The number of ways n tests can sum to s successes, 





Sora! 
N. | ee (2) 
and the value for M., 
M = n+s -2 5 
Ss n-l : ) 
which was just proved, provides the derivation for 6... 
6. _ ds 7 (nts-2)! (n-1l)! s! 
Ng (n-1)! (s-1)! (nts-1)! 
Pp S 
6, = 4 Tec] for s > 0, (4) 
0 for s = 0. 


EUS , 6. is equal to the number of successful cycles, divided by the 


total number of cycles minus one. 


tek k de tek & & & 
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3. Derivation of the Variance. 
The variance of a random variable is defined by 
var (X) = E(X2) - E“(X). 
Therefore, the variance for the unbiased estimator 6. can be 
obtained by solving the following equation: 
var (6.) = E(@,) - E°(8.). 
Equation (4) gives the value for 8. such that E(6.) = 6. Therefore, 
E*(8.) = 94, However, since the derivation of E(6.) provides a 
formula which is useful in deriving the second moment E(6), the 


reverse computation will be shown. 


As - = S i ne 
Ce) Se SEs hy CO) = GS My 0, 
$30 $29 
oe 
- ao $ (y) & 


$20 
Since the s=0 term equals zero, the summation index can be changed, 
and Fors) co 
Ss n+s _ 9° e 9 2 ook 98 
s-l S 
£20 


$29 


By use of the following identity (also derived in Appendix JI), namely 


, 
S(“)e = Gar. 


¢z2O 
a /u-l+s}_s =e 
the summation SG : Je = (i-6)" ; 
$=0 
and 
ail n+s-2\,Ss : n 
SC = GS (6) 
S20 
Thus n 
A (1-0) Q 2 Aa 2 
E(0@,) = ~———— = @, and E(O =i Wels 7 





Similarly, to find the second moment, we note that 


E (6) = > (6.)° N,@ (1-9)" 
$20 
ee 2 
= (1-6)" >(—24} ing a (8) 
$70 











a (=) -(— 





al(aea - + (sea) 
rs Ges ) Eo: 


A ra + (aaoy} 


Substituting this expression into equation (8), and using equation 


It 


(6), we obtain 


ae)" ite *}e - (1-8) > (asa 7 re ¥ 
@- 1-8)" 2 (sie act} (877) e" 


n-1 
] 
Multiplying by a factor gn-1 , we have 


E (64) 








i? ee _ Gay ( —2- nt+s-2 n+s-1 
EE) 2 n+s- at \f n-1 : 
Porm) 
- Se _ -n 1 n+s -2 n+s -1 
5 = ocr) ar SL, (8:2) g 
$?0 
= @- O(n-1) R” £(@) , (9) 
= I 
ie ie 
where {£(9) = y —— eeu e* (10) 





Differentiating equation (10), we have 





044 
1 = = n+s-2 
£(@) = g(@) = D Stet (nts7?) 
$:0 


fo -2) 

n-2 +s -2 

: Zt Jes 
S20 n-l 


Thus, by substituting equation (6), we obtain 


git 
g(0) =  Te)t - 


This expression must be integrated to obtain f£(0®). Thus, 


£(0) 


face) +C = G@)+C, 


where; giv t 


G(®) = Gee” : 


Integrating by parts, letting u =0"°~ and dv = (1-0)7"do, 


n-l fs 2 
n-l \ 1-6 n-2 \ 1-6 n-3 {1-0 


+ eps + (525) + -y"(-8] 4° (-1)) iniGl=e ye a) 


To find the constant C, equations (9) and (11) are evaluated at 


| @ = 0. The constant equals the difference f(@) - G(6) for all @. 


£(0) = 0, and 


G(O) 


+ ine) = 10 








Thus the constant C is zero. ‘Therefore, 


ee re? 2 SR ip bas 
OES n-1 ‘ Hee S “ass : = 5 
+I P SR + (19 R + (-1)" In(l-9). (12) 


Substituting equation (12) into (9) gives the solution for E(@). 


Ao 1 aes: le 
B®) = 6 - em] e-abk + eee 
-1 E : | 
+ (-1) +R a Giye poeas (hears In (1-8) | _ (13) 


Substituting equations (13) and (7) into equation (5) provides the 


derivation of the variance. 


i i ee: 
var (6,) = ofc1-0) = @nty| rR : Ss ec 


3 
— 
j— 


(Cede ; pn-2 4 ce hs CryPeinc-e)]f (14) 


te ke & se te te ke se ok 








4. Numerical Analysis of Variance. 
The equation (14) for Var (6,) lends itself nicely to re- 
cursive numerical analysis. If the alternating series in the 


square brackets is considered, for various values of n, we have 


n=2 R + Rone) = R(1 + Rln(1-8)), 
1 2 3 = 1 ‘ 
n=3 so Re In(1-6) = R¢ ate R(1 + R in(1-8))), 
2 
n=4 =R - =R ay R? Ur RIn(i-@) = 


(15) 
rd - RC F - RU +R In(1-0)))), 


and so forth. If RW, is set equal the sum of the series in equa- 


tions (15), Q is defined as -R, and L = (1 +R 1n(1-0)), then 


: 1 
Cn ae 


I 
, = Qt FZ )Q+R, 


= 
iT 


1 1 1 
We = ((LQ + 5 )Q + 3 Q+ G3 ete. 
Var (6, ) can be re-written as 


var(8,) = @([(1-0) - (n-1) R Wal , 


where Wp is defined by the following recursive formula, 


Wo (1 +R Iln(1-8)), and 


W 
n 


ik 
Ce aan 


10 








Appendix ITI provides tables of var (8.) for 15 values of 9 
and representative values of n between one and 50. ‘input data 
with nine-place accuracy was used in compiling these values. 

Any output error (E,) in this data is a result of the natural 
limitations on input data degrees of significance, and occurs in 


the bracket summation. We have 


al - on 
RE = -==-+4(1)" (6660-246) 2 = (apn Genre 


3 i 


#(-1)""*¢,25) RB" 4 @1ph(.3332--e) ke 


+ (-1)" Rr (inQ-6) + —) . 


Thus, for input data significance of degree 2 3, 


n -3 a6 1 n-J 
EB, = (-1) Re, + R €, +R €, +R 'é, +R Bg-- ete,| 


© 


For values of 02.5 (i.e., for vaiues of R& 1.0), the com- 
puted variance has approximately the same degree of accuracy as the 
input data. For smaller values of 90, however, E, rapidly exceeds 
the value of the variance as n increases. 

As an example, with an input accuracy to nine decimal places, 


Oe= 0.05, R= 19, and ni= 6 


J 


- 6 ae 
E, » (0. €,x 10 J (19 ) = (0.€,x 10 7) (hk, 7R10') 2 (4.7x10 *) . 


but the variance at n = 6 is less than 9 x lor: The tables of 

Appendix II indicate this decrease in accuracy for values of 6<€ 0.5. 
Since @ is normally larger than 0.5 for ttems that are cycled 

to failure, the inherent error term for small @ is expected to place 


little restriction upon the use of the tables. 
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APPENDIX I 
PROOF OF COMBINATORIAL IDENTITIES 
ive To proves: 
nh ° ° 
iy TeX 1 
> cor (fr) = 2 
a2 


X-1 2 
Each ( term can be written as a linear combination of 5") 
n 


terms, where j = 0, l, 2, ---, n, by successive application of the 


(*-*) = ae faeces! 
n n n-1 ; 


rules 


Thus : 


(3) 
(a) > 
a) z fee Siwy 2 es (5=2)("5" 


aN") + (BVP) + (R= == = + (8) 


! 
uum 
3 
| 
= 
ee al 
Fi) 
a 
a 
Oe! 
+ 
nea, 
3 
mI 
— 
a 
a 
eee 
Nod 
20 
meee 
'] 
{ 
8 
i 
+ 
i a 
oS i 
' 3 
i [= 
~~” 
ri 
~ 
3 
rs 
oe 


ke 





5 ta @ | Th 


(to the last term) co? (n\n ys} . 


By summing on like coefficients, we have 


+ FBO arr) 


Az=f | €20 


= 1+ 2(%) 2 1)" \ec| . 


It will now be shown that the summation in the square brackets 


equals zero. 


(a\(r24] ~ i ay eo ey! “ =(i)(z)- 


al) || Geeeal 1 L] |r 
Thus, the summation in the square brackets can be written as 


(? 2 ni () 


EF) 


But this equals zero, as can be seen from 





(x-by)” = (3) x + (7) xr ol yt-e+s-+--e +(7] yn 
Letting x = 1 and y = -l, the solution is provided. 
Ts a r yy ae = i /r 
iL 
(a) -(B)+(5] ~~ + eo" (e) = Zevt (to. 
420 


Je Se iva as 


ibn To prove: 





N+1L}) 1 : N+1) .i 2 I. 
> ( N Je > ( at (i-9) ett. 
tO 20 





This proof can be seen from the successive differentiation of 


2a, 
4-0 
l = 
—— (1 + i) eo 
(1-8) 2. : 
il geo f ‘ 
eee ee 


eee. by induction, 


(1-9) N+ 


m9 


eS SE SS SS Ge Se oe 


14 





ere 


oa) 


o>) 


10 


2 


14 


16 


18 


20 


25 


30 


52 


40 


45 


50 


@ = .999 
00099900 
00000591 


00000099 


less 


than 


107° 


APPENDIX II 


TABLES OF VARIANCE 


6 =92 

00990000 
00036517 
00009262 
00004860 
00003268 
00002459 
00001970 
00001643 
00001410 


00001234 


00000988 
00000824 
-00000706 
00000618 


00000549 


.00000430 
00000353 
00000300 
-00000260 
00000230 


00000206 


@= .95 

04750000 
00558351 
<00E95332 
00109737 
<U0C 756572 
00057524 
00046367 
00038820 
.00085579 


-00029274 


00023491 
.00019614 
00016835 
00014746 


.00013118 


~00010279 
.00008451 
~00007175 
.00006233 
.00005510 


00004937 


15 


@= .90 

09000000 
.01558428 
00653683 
00391053 
-00275400 
00211750 
00171767 
00144401 
00124521 


00109435 


00088071 
.00073676 
00063321 
.00055516 


00049423 


00038780 
00031907 
00027107 
00022207 
00020831 


.00018671 


8 = .80 

. 16000000 
04047190 
01976405 
01258848 
OUD S 7 ley 
00714463 
00585661 
00495849 
00429757 


0087 9S 


00306739 
00257505 
00221867 
00194882 


200173744 


e00136670 
.00112630 
00095780 
.00083315 
-09073720 


.00066197 








n 
LL 


10 


WZ 
14 
16 
18 


20 


25 
30 
SD 
40 
45 


50 


8 = ./0 

- 21000000 
06479650 
03446014 
02284705 
-01694454 
201342257 
01109697 
.00945152 
.00822783 


00728301 


00592068 
00498649 
00430640 
200376929 


~00338292 


00266745 
00220165 
00187430 
00163166 
-00144463 


00129606 


8 = .60 

. 24000000 
08434420 
04754107 
03245893 
02448095 
OL 5g9c! 
01632064 
01397284 
SOE2 Zi, 


01084162 


00885224 
00747794 
00647221 
00570452 


00509941 


00403018 
-00333140 
00283904 
00247343 
sO0ZE91ZS 


-00196682 


@ = .50 

«25000000 
.09657359 
~05685282 
BOSE PAGTIET) 
SLO SLO) HZ oyll 
002453463 
-02055846 
.01768180 
0 l>UGHZ 


sOESSUDL/ 


01131743 
»0095387 26 
00831498 
00734031 


00656988 


,005203'83 
00430779 
00367488 
.00320407 
00284018 


00255049 


16 


e@ = 40 
- 24000000 
09974306 
06077082 
0432656 
0334686 
20272462 
70229504 
AOIGZS 
2017441 


015567 


.012809 
01087 
00945 
00835 


0074 


0059 


004 


Q@ = .30 


. 21000000 


09256908 
.05801098 
.0419615 
O32 7802 
0268734 
2022754 
~OEUT ZS 
017402 


01556 


20025 


009 





10 


rz 


14 


16 


18 


20 


. L6000000 


-07405936 


0475250 
034849 
02746 
02265 
e012 
nOOy, 
014 


~013 


5[0a 


GO = 7.20 


~09000000 
0434201 
.028436 
-02110 
.0167 
BOS. 


SG 


@ = .05 


04750000 
0233439 
015429 
-OLIS2 
0091 


007 


Ke ke ee le tk te ok ke oe oe te 


Ly 


Q@ = .03 


02910000 
~0144022 
200955 


007 


9 = .0l 


.00990000 
-00493 


003 


























